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Multipole electric and magnetic giant resonances in metal clusters are reviewed 
and compared with their counterparts in atomic nuclei. The main attention is paid 
to El resonance (dipole plasmon). 

1 Introduction 

Metal cluster (MC) is a bound system consisting of atoms of some metal. 
The amount of atoms can vary from a few to many thousands. Some MC, 
mainly of alkali (Li, K, Na, ...) and noble (Ag, Au, ...) metals, demonstrate 
a striking similarity to atomic nuclei (see reviews^"^). In these clusters the 
valence electrons are weakly coupled to the ions and, like nucleons in nuclei, are 
not strongly localized. The mean free path of valence electrons is of the same 
order of magnitude as the size of the cluster. This favors the valence electrons 
to form a mean field of the same kind as in nuclei (with the similar shell 
structure and magic numbers). In addition to the mean field, MC demonstrate 
other similarities with atomic nuclei: deformation in the case of open shells, 
variety of giant resonances (GR), fission, etc.. As a result, many theoretical 
ideas and metkojis of nuclear physics can, after a certain modification, be 
applied to MCHdM 

This review is devoted to collective oscillations of valence electrons in MC. 
Valence electrons can be considered as the counterparts of nucleons in nuclei, 
and their oscillations as the counterparts of nuclear GR. Investigation of GR 
in MC is interesting in two aspects: it allows to understand deeper general 
properties of collective modes in finite Fermi systems and, simultaneously, al- 
lows to study striking peculiarities of MC. GR in clusters and atomic nuclei 
are well overlapped. However, some specific properties of MC cause consider- 
able differences in the behavior of GR in these two systems. For example: the 
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Coulomb interaction and the " spill-out" effect provide a specific dependence of 
GR properties on the mass number; the negligible character of the spin-orbital 
interaction leads to the decoupling of spin and orbital magnetic modes; clusters 
can have much more particles (atoms) than nuclei, which favors very strong 
orbital magnetic resonances; for most of the clusters the role of the ionic sub- 
system is important; at different temperatures MC can be in solid, liquid and 
even "boiling" phases, which greatly influences GR properties; characteristics 
of GR vary considerably whether the clusters are charged or neutral, free or 
embedded to a substrate, pure or with impurities atoms, etc.. 

Our consideration will be limited by certain physical conditions. 

- i) The modern techniques allow to fabricate atomic clusters from atoms of 
about any element of the periodic table. However, the conception of the mean 
field for valence electrons is realized only for a minority, - mainly for clusters 
of alkali and noble metals and, in a less extent, for neighboring elements. So, 
we should limit ourselves by this MC region. 

- ii) In some alkali metals (Na and K) the ionic lattice can, to good accu- 
racy, be replaced by a uniform distribution of the positive charge over cluster's 
volume. This is so-called jellium approximation which greatly simplifies the 
analysis and calculations. This approximation is enough for the description of 
many properties of alkali MC and will widely be used in the review. However, 
it often fails beyond Nsuand K and then a more explicit treatment of the ionic 
structure is necessary 

- iii) The ionic subsystem is supposed to be " frozen" . In the review only 
collective oscillations of valence electrons will be considered. 

- iv) The validity of jellium approximation is supported by temperature 
fluctuations of ions, which smooth ion positions. It fails in the low temperature 
region (approximately at T < 100 K) where the explicit treatment of the ionic 
structure is important. At too high temperatures {T > 1000 K) the quantum 
shells of the mean field are washed out, what establishes an upper limit for our 
considerations. We will consider GR in a temperature interval between these 
extreme cases. 

2 Theoretical Grounds 

Due to the similarity between MC^nd nuclei, many models of nuclear theory 
have been applied to study MCd'Q. In due time, some of them have been 
introduced to nuclear physics from solid body field and then subsequently 
modified to describe finite Fermi systems. Now they turn out to be useful for 
clusters. In particular, a large variety of the RPA methods have been adopted , 
scaling from simple versions, like the sum rule approach^"^ and the local RPA 
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00, to sophisticated full RPA models, like time-dependent Hartree-Fock 
and time-dependent local densityj-ajiproximation (TD-LDA)^^~^^ (for a more 
complete list of citation see RefsoaO). The simple models can describe the 
gross structure of GR but not the fragmentation of the collective strength. The 
full RPA models can describe the fragmentation but are very time consuming. 
The last shortcoming becomes crucial for deformed and large spherical clusters 
where the number of particles, and thus the size of the configuration space, is 
very large. In this connection, the intermediate class of the models, the RPA 
with separable residual forces (SRPA), seems to be very promising^'^^"^^. The 
separable ansatz allows one to turn the RPA matrix into a simple dispersion 
relation. This drastically simplifies the eigenvalue problem preserving, at the 
same time, the main advantage of the full RPA to describe the fragmentation 
of the collective strength. The SRPA iiersion derived in Refs.^^"^'^ provides 
the accucajCjf|-Qf full RPA calculations cil, can be applied to systems of any 
shapeEa'E3E3S, and-aLLoats to treat GR in MC and atomic nuclei on the same 
microscopic footingc3E3c3. The results obtained within this SRPA version will 
be widely used in the review as illustrative examples. 

For the description of collective oscillations, the SRPA, and most of the 
u^r models, exploit, as a starting point, the Kohn-Sham energy functional 
for a system of N^. valence electrons: 

E{n{r,t),m{r,t),T{r,t)} = 1/2 J T{r,t)dr + J w^c(ri(r, i), m(r, i))dr 

+ 1/2 J J ("(r.0""4r))Kri,0-n.(ri)) ^^^^^^ 

which includes the kinetic energy_th£ exchange-correlation term in the lo- 
cal density approximation (LDA) and the Coulomb interaction, respec- 
tively. Here, n{r,t) = n{r,t) 1 +n{r,t) J,= Si l<^i(r, m{r,t) = n{r,t) 1 
—n{r, t) I and T(r, t) — J2i I V 4'i{^: OP are the density, magnetization density 
(z-component) and kinetic energy density of valence electrons, respectively; 
ni(r) is the ionic density in the jellium approximation; (j)i{v,t) is a single- 
particle wave function. The convention e = me = ?i = 1 is used. The func- 
tional (1) can have additional terms if the ionic structure is treated beyond 
the jellium approximation. 

The time-dependent single-particle Hamiltonian is obtained as 

H{v,t)Uv.t) = -^ (2) 
64>i (r,t) 

In the small-amplitude limit of a collective motion, the densities can be written 
as n{r,t) = no{r) + 6n{r,t) and m(r,t) = mo(r) + 6m{r,t) where no(r) and 
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r7io(r) are the static ground state densities (mo(r) = in spherical and unpo- 
larized clusters) and the values Sn(r,t) and (5m(r,i) are small time-dependent 
density variations (transition densities). Then, in the linear approximation to 
the density variations, the Hamiltonian (2) is a sum of the static and dynamical 
parts. The static part 

Ho{r) = T + Voir) = -y + (-^)„=„o,m=mo + J dr^ (3) 

constitutes the Kohn-Sham single-particle potential. It can be approximated 
with a good accuracy by phenomenological potentials_siich as the harmonic 
oscillatofi (for snaajlLpilierical MC), Nillson-Clemengeiocil (for deformed MC) 
or Woods-SaxonEifE3c3 (for spherical and deformed MC . 

In the electric channel, the dynamic part of the Hamiltonian (residual 
interaction) has the form 

5Hir,t) = (^?!^)_„„5n(r,i) + / P^dr,. (4) 



dn^ J |r — ri 

The dominant term here is the Coulomb interaction. The residual interaction 
in this channel is always positive (repulsive) and shifts the unperturbed electri- 
cal multipole strength from the typical particle-hole (ph) values ujph — 0.9— 1.5 
eV to higher energies 2.6-3.2 eV. 

In the spin channel, the dynamical part is 

dH{r,t) = ( )n=no,ni=moSn{r,t)dm{r,t). (5) 

anam 

Here, the residual interaction is determined by the exchange-correlation term. 
It is always negative (attractive) and shifts the unperturbed magnetic multipole 
strength from coph = 0.9 — 1.5 eV to lower energies 0.2-0.8 eV. The residual 
interactions (4) and (5) are typical for the TD-LDA calculation scheme. 

In what follows we will mainly consider clusters constituted from monova- 
lent atoms, like alkali metals, for which the numbers of valence electrons and 
atoms coincide, Ne = N. 



3 Electric Dipole Giant Resonance (El GR) 

Unlike nuclei, where different kinds of GR are well investigated both exper- 
imentally and theoretically, our knowledge in clusters is mainly limited by 
the electric dipole resonance (dipole plasmon). Experimentally the El GR has 
been observed in a variety of clusters: small and large, spherical and deformed. 
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neutral and charged, hot and cooled (see, for example, Refs.^^^'*^). As a rule, 

the photoabsorption cross section was measured by methods of the depletion 

spectroscopy. For other GR {EL{L ^ 1), ML) there are only theoretical 
predictions6-9'28.30,32,50.5i^ 

Physical interpretations of El GR in clusters and nuclei are very similar: 
while in nuclei it is caused by translations of neutrons against and protons, 
en in clusters it is a result of translations of the valence electron against ions 
. In spite of this similarity, the dipole resonance in clusters exhibits many 
interesting peculiarities which will be discussed below. 

3.1 Energy of El GR: Step by Step 

The description of El energy for clusters is a rather complicated task. For 
example, while in nuclei its energy depends on the mass number as A~^^^, 
in clusters the El-energy can both decrease (Ag clusters) and increase (alkali 
MC) with the number of atoms. Let us consider this important characteristic 
step by step. 

Step one: Mie frequency and spill-out effect. In the simplest ap- 
proximation, MC can be considered as a^assical metallic drop. Then, the 
El-energy is described by Mie expression LJ: to Mie = Wp/v^ where Up is the 
plasma frequency. For Na clusters iOMie = ?)AleV . This value is much higher 
than the experimental El-energy which is 2.5-2.8 eV for spherical Na clusters 
with N < 100. 

The agreement with the experiment is considerably improved if we take 
into account the quantum spill-out effect. This effect means that since the 
valence electrons are quantum entities, they are not well localized and so, 
unlike the classical ionic jellium, can be partly spilled out beyond the jellium 
boundary. In principle, this effect takes place in any two-component quantum 
system including atomic nuclei and atoms (a "neutron skin" in small nuclei is 
a relevant example). With the spill-out, the El energy in MC is described as@ 

^^El=^^Ane(l-^—) (6) 

where SNe is the number of spilled out valence electrons. As a result, the 
discrepancy with the experiment reduces to 0.2-0.3 eV. The spill-out effect 
allows to explain the increase of the El-energy with N, observed in alkali 
MC. The value 6Ne decreases with the size_(for example, 6Ne = 1.5(19%) 
and 9.5(7%)in Nag and Nai^s, respectively o) leading to the corresponding 
increase in the El-energy. 

Step two: beyond jellium approximation, ionic structure, local 
and nonlocal effects. The remaining discrepancy can be removed in a large 
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extent by the explicit treatment of the ionic subsystem. First of all, we should 
take into account that ions are not the points but have a size. Inside this 
size, ion core electrons (ICE) (do not confuse them with the valence electrons) 
screen the pure Coulomb interaction of ions and valence electrons. To take into 
acroiint this scrppning the atomic pscudopotcntials (PP) are used (see, for in- 
stance, CJ'ESE^Iij). They allow to describe correctly the spectrum of valence 
electrons in isolated atoms without the solution of the complicated many-body 
atomic task. Being a sum of contributions of ICE with different orbital inn- 
menta, PP have local (s-electrons) and nonlocal (p and d electrons) parts 
To avoid dealing with nonlocal functions^ |tie Pseudo-Hamiltonians (PH) were 
introduced as the next simplifying stepEaH. PH, being derived from PP, lead 
to less involved (but with the same accuracy) calculations since, unlike the 
PP, they treat the nonlocality only through the differential operators. Fold- 
ing atomic PH with jellium, one gets PH for atomic cluster s^'^ . PH have 
the additional advantage to be easily incorporated to the common calculation 
schemes. 

As compared to the conventional Kohn-Sham Hamiltonian, PH include the 
additional local and two non-local (the orbital contribution and the effective 
mass) terms. As is seen from Figure 1, in K clusters (the same for Na) the 
nonlocal contributions are negligible and the local term is enough to get good 
description of the El-energjcJ. This is not the case for Li clusters, where 
only both, local and nonlocal, contributions proszide the agreement with the 
experiment EJ. In some studies (see, e.g., Ref. E3) the ICE effects are taken 
into account together with some averaged treatment of the ionic arrays in a 
cluster. The latter leads to an additional, but rather moderate, redshift of the 
El-energy. 

Step three: direct dynamical ICE contribution. The ICE effects 
discussed above are realized through the change of the single-particle charac- 
teristics with the subsequent renormalization of the residual interaction. Be- 
sides this indirect way, the ICE can directly influence the dynamics and thus 
lead to new peculiarities of the El GR. This can be well demonstrated for 
Ag clusters where, dike atomic nuclei and unlike alkali MC, the El-energy de- 
creases with a sizeCJ. The physics behind is that in these clusters the energy 
of ICE excitations is comparable to the El-energy. The coupling of these two 
modes additionally screens the interaction of valence electrons with ions (di- 
rect dynamicaJ-JCE contribution) and finally causes the redshift (decrease) of 
the El-energy e2I. This effect is mainly of a volume character and so intensifies 
with a cluster size. As a result, the El energy in Ag clusters decreases with N. 
This tendency overpowers the opposite one caused by the spill-out effect. 

In alkali clusters, where the ICE excitations have energies much higher 
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than the El GR, the direct dynamical ICE contribution can be neglected and 
the evolution of the El-energy with the cluster size is determined mainly by 
the spill-out effect. 

4 Landau Damping and Width of El GR 

The main physical mechanisms forming the plasmon width are the thermal 
fluctuations of a cluster shpic_and^iia Landau damping (RPA fragmentation 
of the collective strength) qEjoOl^O'EJ. The relative contributions of these 
two mechanisms change with a cluster size. As is seen from Figure 2, in 
small clusters, like Nd^i, where the Landau damping is negligible, the thermal 
fluctuations determine about all the width. In clusters of moderate size, like 
iVor^g, the Landau damping is stronger and greatly contributes to the width. 
This is especially the case for deformed clusters. In large clusters, like Na'^^-^, 
the Landau damping is weaker though its contribution to the width remains 
to be considerable. 

Figure 2 (bottpm) shows that the Landau damping is closely related with 
the shell structured. In Na^i the dipole plasmon lies in the wide gap between 
the bunches of ZW =1 and ZW =3 particle- hole (ph) states and remains almost 
unperturbed as a collective peak. With increasing the cluster size, the reso- 
nance approaches the bunch AA/'=3 and, in A^a^g, is already interferes with 
ph states of this bunch, which leads to the considerable Landau damping. For 
larger clusters, the plasmon runs to the swamp of ph states. This leads to a 
general trendy of. increasing the width which is, however, overlaid by sizeable 
fluctuationsEja. But here a further mechanism comes into play: the coupling 
between the resonance and ph states fades away due to increasing mismatch of 
AA/'=1 ph configurations (which mainly generate the plasmon) and surround- 
ing ph states with much larger values of AA/". This finally leads to a decrease 
of the plasmon width cx A^e ^^'^ estimated analytically in the wall formula 
and tested in the RPA calculationsE3. 

The Landau damping in MC with A^ < 40 is rather sensitive to clus- 
ter charge: being strongest in negatively charged ones (anions), the Landau 
damping is considerably rediiced while passing to neutral and then to positively 
charged clusters (cations) E3. This effect is caused by a strong dependence of 
the single-particle potential depth Vo on the cluster charge. In anions the po- 
tential is shallow (Vb — —2 eV), the energy gaps between AM bunches are 
very smooth and, so, there are good conditions for a sizeable Landau damping 
(see discussion above). In neutral clusters and more in cations, the potential 
depth is increased to about -7 eV, the gaps between AA/" bunches in the ph 
spectrum become more distinctive, which weakens the Landau damping. 
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5 Temperature Effects 



In most of experiments with GR in MC, the typical cluster temperature is 
estimated to be in the interval 300-900 K which corresponds to the thermal 
energy kT = 0.03 — 0.09 eV. At these temperatures, ions behave as classical 
particles and quantum properties of the cluster are mainly determined by va- 
lence electrons. This can be easily provedE3 by using the uncertainty relation 
AxAp > h. This relation gives lower bounds for the momentum and energy of 
a particle in a system: Ap = h/ Ax and AE = (Ap)"^ /2m, respectively. Taking 
Aa; < l.bA (the diameter of Na2()) for both ions and valence electrons, one 
gets 

AEe > O.WeV for electrons, 
AE.^ > IQ-'^eV for ions. 

The energy of a quantum motion of valence electrons considerably exceeds the 
thermal energy, which favors their quantum behavior. The opposite situation 
takes place for ions which, therefore, should exhibit the classical behavior. 
Such result takes place because the ionic mass is much larger than the electron 
one. 

The difference in ionic and electron masses leads to other interesting con- 
sequence. Namely, almost all the thermal energy is contained in the ionic 
subsystem. Valence electrons are embedded to the thermal ionic bath. So, 
unlike atomic nuclei, MC represent the case of the canonical ensemble. 

The bulk melting points for K, Na and Li are Tb = 336, 371 and 452 K, 
respectively. This means that most of the measurements for GR in MC have 
been done for clusters in a liquid-like phase. 

As was mentioned above, in small clusters, thermal shape fluctuations 
provide the dominate contribution to the plasmon width. It is interesting 
that, while in nuclei these fluctuations are mainly of a quadrupole form, in 
MC they are mainly octupoleO. The reason is that spherical and neighboring 
MC are soft to the octupole deformation. 

Photoexcitation is a rapid process in the ionic time scale. So, every re- 
sponse of a cluster represents its instantaneous shape and the expepimental 
cross section gives a properly weighted response of all allowed shapes l3. 

The higher the temperature, the larger the plasmon width and the smaller 
the plasmon energy. The lenmerature shift is estimated as about 1% of the 
plasmon energy per 100 Kou. It can be explained by the effectively increase 
of the cluster size with temperature. The larger the size, the bigger the static 
dipole polarizability which is expressed through the cluster radius as a^i = R^. 
The polarizability is connected with the plasmon energy through the inverse 
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sum rule, aEi = 2rn_i ~< El >^ /(^ei- So, the higher the temperature, the 
larger asi and, consequently, the smaller ujei- 

Recent experiments show that at sufficienily low temperatures the gross- 
structure of the El GR drastically changesllj. For example, the axially de- 
formed cluster -/Va^*! at 380 K demonstrates the typical two-peak spectrum 
determined by the deformation splitting of El GR. At 35 K the same res- 
onance exhibits much more complicated structure including at least 6 well- 
distinguished peaks. This structure reflects the ionic arrangement which is 
not washed out at so low temperature by variations of ions. In this case, the 
jellium approximation is not valid and models based on this cannot be applied. 
The El GR in small clusters at low tempacature seems to be best described 
by ab initio quantum-chemical calculations eJ. 



6 El GR in Deformed Clusters 

Like nuclei, MC with open shells have quadrupole deformation'*^"'*^'^'*"^^. 
There are experimental indications of both prolate and oblate axial quadrupole 
shapes, as well as of 7-deformation^^~'*^. In the framework of different meth- 
ods (Strutinski's shell correction method, ultimate jellium model, etc.) hex- 
adecapole and octupole deformations as well as high isomerism have been 
predicted^''"^^. Rather strong quadrupole, hexadecapole and. octupole defor- 
mations should take place at least up to MC with N ^ 700E3. Like in nuclei, 
El GR in axially deformed MC exhibits the deformation splitting in two peaks 
(see Figure 3). The right peak is about twice larger than the left one in prolate 
clusters (see Nafi, Naf^, Na^-j) and, vice versa, in oblate clusters (see Nd^^). 

Most of MC are deformed. But getting an experimental information on 
a cluster shape, even in the simplest case of a quadrupole deformation, is 
rather nontrivial problem. In nuclei rotational bands serve as source of such 
information. In principle, deformed clusters can rotate. But, due to a large 
value of the moment of inertia, rotational energies are very small and, being 
of the same order of magnitude as the thermal energy, fail to be observed. In 
this connection, the splitting of El GR in deformed clusters is now a single 
direct manifestation of quadrupole deformation and the main source of the 
information about it. 



7 Multipole GR, Asymptotic Trends, Restoring Forces 

So far, the depletion spectroscopy methods (photoabsorption and photofrag- 
mentation) were mainly exploited for observation of El GR in MCB. The other 
reactions ((e, e'), (7,7') and etc.) are not yet sufhciently developed, which im- 
pedes the observation of other GR. The similar situation took place in nuclear 
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Table 1: Relative contributions to 7713 for Na%2'- kinetic energy (7713 (T)), exchange and cor- 
relations (m3(a;c)), electron-electron interaction (mz{ee), eJectron-ion interaction (ms{e{)) 
and total Coulomb interaction (7713(C) = 7773(66) + 7773 (67)) Q. 



L 


7773 (T) 


''773 (xc) 


7773(66) 


7773(67) 


7773(C) 


1 











1 


1 


2 


0.08 





-0.77 


1.69 


0.92 


5 


0.51 





-2.29 


2.78 


0.49 



physics in early seventies. For this reason an investigation of EL GR with 
L 7^ 1 is yet hmited to theoretical predictions^"^^^*'^"^^^'^"'^-^. In Figure 4, E2 
and E3 GR in spherical Na^^, calculated within the SRPA, are presented as 
typical examples. 

It is instructive to consider the main trends of electrical multipole giant 
resonance with the size (N) and multipolarity (L), and also the origin of the 
GR restoring forces. Such analysis has been done within the sum rule approach 
(SRA) in Refill. In the jellium approximation for valence electrons, no(^) ~ 
ni{r) = n^6{r — R) (tbe spill-out effect is neglected), the energy of EL{L ^ 0) 
GR can be written asQ 



where mi = B{EL, gr —^ i)uJi and = B{EL, gr —> i)ujf are the sum 
rules, (3f = (3/5)^/2(37^2)1/3^^1/3 ^ ^j^^ radius of a cluster. The first 
term in Eq. ^ is the contribution of thje kinetic energy (the similar expression 
have been obtained earlier in Ref. O). The second term is determined by 
the Coulomb interaction between valence electrons (ee) and valence electrons 
and ions (ei). Eq. |^ shows that El GR is determined only by the Coulomb 
interaction. In the limit of large R, one has 



The larger L, the higher the excitation energy of the GR. In general, due to 
the first term in Eq. |^, the energy of EL{L ^ 0, 1) GR is decreased with N. 
For low L in small clusters this tendency is changed by the spill-out effect. 

The separate analysis for EO GR gives the increase of the EO energy with 
N to the limit ojeo — > fiUn. 

IL-is seen from Eq. R that the value has the meaning of a restoring 
forced. In Table |l| the contributions to TO3 from different terms of the Kohn- 



10 



Sham functional (1) are presented. It is seen that the restoring force for El 
GR is determined by the electron-ion interaction (ei) only. With increasing L, 
the electron-electron contribution (ee) raises and starts to compensate the (ei) 
Coulomb part. Simultaneously, the kinetic energy term grows. For high L, the 
contribution of the total Coulomb interaction goes to zero and all the restoring 
force is determined by the kinetic energy. The purely volume exchange corre- 
lation term (xc) which within the LDA depends only on the electron density 
does not contribute to 7713. 

The restoring force should not be confused with the residual interaction. 
As is seen from Eq. ^, the residual interaction, unlike the restoring force, has 
only the (ee)- and (xc)-terms for any L (where the (ee)-term dominates). 

8 Anharmonicity and Multiphonon GR 

How much harmonic are the GR in metal clusters? How strong is the mix- 
ing of one and two phonon states? Investigations performed within different 
approaches^^"^^ give contradictory answers for one-phonon GR. While the 
shCjlLmodel calculations found for El GR in A^a2o some signals of anharmonic- 
ity Il3, Oth|£i; studies predict the harmonic character for M2(spin-dipole) and 
EL GRIIjO. It should be noted that these studies have been done for rather 
small clusters with iV < 20. In this size region the GR energy lies safely below 
the lowest 2p-2h configurations, what does not favor the anharmonic effects. 
This picture can change in larger clusters where GR approach the region of 
2p-2h configurations.!—. 

The calculations O predict a noticeable anharmonicity for most of double 
(two-phonon) GR placed at 8-15 eV. These GR exhibit a weak mixing with 
one-phonon states but a considerable fragmentation between two-phonon con- 
figurations. Most strong effect is expected for some 0+ double GR, for example, 
for (1~ (8) l~)o+ in ^O'ti- With the appearance of new experimental techniques 
allowing investigation of multiple GR these predictions are quite important. 
The techniques use non-intenseJemtosecond laserslla or exploit collisions of a 
cluster with highly charged ionsES. Quite recently the-caultiple GR constructed 
from 3-4 dipole plasmons has been observed in iVa^Jgllj. 

9 Magnetic GR 

Like in atoms, the spin-orbital interaction in metal clusters is negligible and 
thus spin and orbital collective magnetic modes are well decoupled. The sep- 
aration of these two modes in MC is easier than in nuclei. 
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Table 2: Relative contributions to for Nag2 and Nag\2 : kinetic energy {m^{T)), 
correlation {rn^(c)X^a.nd total Coulomb interaction (013 (C)). The data are extracted from 
the Fig. 1 of Ref. 
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Na<32 






m3(c) 






m3(c) 




1 





0.77 


0.23 





0.30 


0.70 


2 


0.50 


0.49 


0.01 


0.31 


0.68 


0.01 


5 


0.76 


0.20 


0.04 


0.59 


0.38 


0.03 



9.1 Spin-Multipole GR 

Magnetic multipole resonances (ML) of spin character caused by theexternal 
field Ql = r!jYLo<Jj were studied within the SRA and RPAi'B0li. For 

L = 1 the operator Qi ~ ^j'^j provides the opposite shifts of spin-up and 

spin-down electrons in z-dircction. Unlike the electric resonances , the residual 
interaction for ML GR is defined only by the exchange and correlations ((xc)- 
term) since only the (xc)-terni depends on the magnetization density (see Eq. 

. Therefore the study of ML resonances can provide a valuable information 
about (xc)-effects in clusters. 

Approximating the electron density by the—expression no = noo/(l + 
exp{{r — R)/a), one gets the energy for ML GRIIj 

+ -i(»r'("o. ».o) - vi«>(-a„,n,„))^fl^ (9) 

m ban 

where v^Jc'\no,mo) = ^^wic'"' (ri, m) |(„=„o,m=mo)- For other notation see 
Eq. |. As compared to Eq. | for EL GR, Eq. (9), was derived taking into 
account the spill-out effect. Furthermore, due to the presence of the spin in 
the operator Ql, the (xc)-term now contributes to m^, unlike the case of EL 
GR. However, the exchange contributions (Pauli principle) to and v^l are 
the same and then, only correlation effects enter Eq. ^. The energies of spin 
ML resonances decrease with N and go to zero for large sizes. The larger L, 
the higher the GR energy. The behavior of ML GR much depends on the 
diffuseness parameter a. 

Table || demonstrates that the restoring force for spin- multipole GR- differs 
from the electrical GR case. Namely, the contribution of correlations^ domi- 
nates for L = 1 and 2 and remains to be considerable for larger L. The correla- 
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Table 3: The excitation enen gii of id strength (within the interval 0-1 eV) of orbital Ml GR, 
calculated within the SRPAe3cj. See the text for notation. 
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n term includes long-range RPA correlations short-range correlations 
and others. The correlations greatly influence both static and dynamical 
characteristics of MC^^'^^'''^"^^ and their investigation is very important. 

9.2 Orbital GR 

Since the number of atoms in MC can be much more than the number of 
nucleons in nuclei, much larger values of the single-particle orbital moment 
can be achieved what can give origin to very strong orbital magnetic multipole 
resonances, orbital ML GR. That is, clusters, as nuclei, can exhibit orbital 
ML GR Jike, "scissors" , twist mode, etc., nevertheless with a much stronger 
strength E. 

Investigations of the specific low-energy orbital Ml GR, which can ex- 
ist only in deformed clusters demonstrated that it can serve as a good indi- 
cator of the cluster quadrupole deformation 7, 81-83 ^ Indeed, in some cases 
the deformation splitting of El GR is washed out by other effects and is not 
enough distinctive to get a reliable information on cluster deformation. Then 
the orbital Ml GR can be used for this aim. Macroscopically, this resonance 
is treated as small-angle rigid rotations of the ellipsoid of valence electrons 
against the ionic ellipsoid. Such collective mode was shown ta^ba coupled with 
the quadrupole component, v(yz), of the displacement field&Ej. The orbital 
Ml GR has the counterpart in deformed nuclei, well known as the "scissors" 
mod(E3. The latter describes the rotations of the neutron ellipsoid against the 
proton one. The orbital Ml GR is represented by = 1+ states {K is the 
angular-momentum projection) with a low excitation energy and strong Ml 
transitions to the ground state. For Na clusters these characteristics are esti- 
mated asiS OJMI = 4.6/32 A^e~^^^(l + 5^)~^/^ eV and B{M1) = Ll/^jiVe^Vh 
where /?2 is the deformation parameter, B{M1) is the reduced transition prob- 
ability and ujo is the harmonic oscillator frequency. Both lumi and B{Ml) are 
proportional to the deformation parameter and so the orbital Ml GR survives 
only in deformed clusters. 
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The results of the first reaUstic RPA calculations for orbital Ml GrE10 
are given in Table 3. It is seen that this resonance has low excitation energies. 
The most remarkable result is that already in clusters with about 300 atoms, 
the orbital Ml GR strength reaches very high values, 700-800 nl- This GR is 
described in detail in RefEj of the present Proceedings. 

10 Other GR in Atomic Clusters 

As compared to nuclei, atomic clusters provide many specific manifestations 
of El GR. For example, clusters embedded in a dielectric matrix demonstrate 
a strong screening effect: the matrix screens the residual interaction between 
valence-electrons in a cluster, which results in the considerable decrease of El- 
energylSa. In mixed and coated clusters the impurity (or coated) atoms much 
influence both the ground state and properties of El GR (see, e.g. Refs.®^~^*). 
In the fuUerene Cgo, two El GR are known as determined_bcy weakly bonded 
TT electrons and strongly bonded a electrons (see, e.g. Ref. Ej). 

'^He and ^He clusters representing collections of fermions {'^He atoms) 
and bosons {'^He atoms), respectively, should be mentioned. In ^iJe clusters 
iasXJLHe atoms (not valence electrons) form a mean field with quantum shells 
These clusters are characterized by strong surface effects. Unlike nuclei 
and MC, ^He clusters represent the case of one- component Fermi-system and, 
so, have no El GR. At the same time, the study of other EL GR reveals new 
possibilitiiiST-for instance, the comparison of the GR properties in Fermi (^He 
clusters) Ej'O and Bose {^He clusters^ systems. 

Summary 

Giant resonances in atomic clusters have been observed. Being much similar 
to their counterparts in atomic nuclei, GR in MC demonstrate, at the same 
time, numerous exciting peculiarities. The unique situation takes place now 
in many-body physics where, in addition to atoms and atomic nuclei, a new 
family of small Fermi systems (MC, fuUerenes, He^ clusters, quantum dots) 
appears. This greatly enlarges our possibilities in many-body studies. All 
mentioned systems possess, in a different extent, a mean field with quantum 
shells. 

It should be noted that atomic clusters are attractive both for fundamen- 
tal studies and practical applications Ej. Last achievements (creation of new 
materials, machinning superhard surfaces, creation of extremely large energy 
densities in a matter, catalysis, microelectronics, microcomputering, etc.) show 
that, due to atomic clusters, one may expect in a recent future a remarkable 
progress in many high-tech fields. 
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FIGURE CAPTIONS. 

Figure 1. El GR (dipole plasmon) in and Lz^j^ calculated in the 
framewori_of the SRPA with (down) and without (up) the nonlocal ICE con- 
tributiona. For Li^j^ the photoabsorption experimental dataE3 (A) in /N^ 
are compared. 

Figure 2. Toul El GR in spherical Na clusters from different size regions. 
The SRPA resultsEil are shown as bars for every RPA state to demonstrate the 
Landau damping and as smoothed by a Lorentz function (of the width 0.25 
eV) to simulate the typical thermal broadening of the plasmon. The length 
of the bars is rescaled by the factor 1/2.55 to fit the scale of the smoothetl 
strength. The photoabsorption experimental data are taken from Ref. E3. 
Bottom: The number of dipole particle- hole configurations, as a function of the 
energy, corresponding to AA/" =1 (light dotted bricks), AA/'=3 (dashed bricks), 
AA/'=5 (dark dotted bricks) and AA/" >7 (unfilled bricks) dipole transitions in 
Na clusters presented in the top of the figure. M is the principal shell quantum 
number. The arrows mark centroid energies of the plasmon. _ 

Figure 3. El GR in deformed Na clusters. The SRPA resultsE3 (curves 
and bars) are given_by the same way as in Figure 2. The experimental data 
are taken fmm RefE2l. The deformation parameter j32 is extracted from the 
experiment 113 following the prescriptiom. 

Figure 4. E2 and E3 GR in Aa+ calculated within the SRPA. 
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